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Several analytic approaches predict for SU {Nc) Yang-Mills theories in Landau gauge an enhanced 
ghost propagator G{p^) and a suppressed gluon propagator D{p^) at small momenta. This pre- 
diction applies to two, three and four space-time dimensions. Moreover, the gluon propagator is 
predicted to be null at p = 0. Numerical studies by several groups indeed support an enhanced 
ghost propagator when compared to the tree-level behavior 1 / and a finite infrared gluon prop- 
agator. However, the agreement between analytic and numerical studies is only at the qualitative 
level in three and in four dimensions. In particular, the infrared exponent of the ghost propagator 
seems to be smaller than the one predicted analytically and the gluon propagator seems to dis- 
play a (finite) nonzero value at zero momentum. It has been argued that this discrepancy might 
go away once simulations are done on much larger lattice sizes than the ones used up to now. 
Here we present data in three and four space-time dimensions using huge lattices in the scaling 
region, i.e. up to 320^ at j3 = 3.0 and up to 128^ at j3 = 2.2, corresponding to V « (85 fm)-' and 
V « (27 fm)^. Simulations have been done on the IBM supercomputer at the University of Sao 
Paulo. 
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1. Introduction 



In the Gribov-Zwanziger confinement scenario for Landau gauge ^ tiie gluon propagator 
D{p^) is predicted to be infrared (IR) suppressed at small momenta. In particular, one should have 
D(0) = 0, implying that reflection positivity is maximally violated. This violation of reflection 
positivity may be viewed as an indication of gluon confinement At the same time, the Gribov- 
Zwanziger im ^ and Kugo-Ojima |^ confinement scenarios predict (in Landau gauge) a ghost 
propagator G{p^) enhanced in the IR limit. This represents a long-range effect and could be related 
to quark confinement [[T], |3[ |[ |||. 

Several analytic studies [0, ||, ^, 10 1 agree with the above scenarios predicting, for small mo- 
menta, a gluon propagator D{p^) oc ^j^^ ^ ghost propagator G{p^) oc 1/ ^2(i+aG)^ ^jth the 
relation ajj = lac + (4 — t/) /2. Here d is the space-time dimension. Clearly, if > 1 one has 
D(0) = 0. In the four-dimensional case, one finds [||, ^ ~ 0.59 and ao = lac- Similar power 
behaviors have also been obtained for the various vertex functions of SU (Nc) Yang-Mills theories 



11, 12] 



These results have been confirmed at the quantitative level in the two-dimensional case, using 
lattices up to almost (43 fm)^ [|I3|]. In the three-dimensional case [14], one clearly sees an IR- 
suppressed gluon propagator. However, using lattice volumes up to about (24fm)^, it was not 
possible to control the extrapolation to infinite volume. In particular, the data for the rescaled 
gluon propagator at zero momentum D{0) could be fitted [as a function of the inverse lattice side 
I/L = 1 /{Na)] using the Ansatz d + b/LF both with d = and with J / [ITl]. Here a is the 



lattice spacing and N is the number of lattice points per direction. Finally, in four dimensions, the 
gluon propagator is clearly less divergent than in the tree-level case [|l^, [l^, 17, 18]. On the other 
hand, even using lattices with a lattice side of about lOfm, one does not see a gluon propagator 



decreasing at small momenta []18[]. One should stress, however, that the Landau gluon propagator 
clearly violates reflection positivity, in two, three and four space-time dimensions [13, 19, EO]. 



For the ghost propagator, the IR exponent ac obtained using analytic studies has been confirmed 
in 2d [13], while in the 3d [ pT] ] and in the 4d [16, 17, 22, 23] cases the exponent obtained using 
lattice numerical simulations is always smaller than the one predicted analytically. Let us also 
recall that, in Ref. [p4{], it was shown that gluon and ghost propagator for SU (2) and SU (3) Yang- 
Mills theories are in very good agreement from momentum p \ GeV to about p k, \0 GeV. 



Similar results have been presented in [25]. These findings suggest that the IR behavior of these 



propagators is independent of the gauge group SU {Nc), as predicted analytically [0, [sj, ^. 



Finally, by solving Dyson-Schwinger equations on a finite four-dimensional torus [12, 26, 27] 
one can show that the gluon propagator (at small momenta) seems to diverge for volumes up to 
about (8fm)^, develops a plateau for V ^ (9fm)^ and is IR suppressed for a lattice side larger 
than 10 fm. Also, the extrapolation of these results to infinite volume gives a null gluon propagator 
at zero momentum. At the same time, one obtains that, after eliminating (for each volume) the 
data corresponding to the first two non-zero momenta, the ghost propagator shows a power-law 
behavior and, in the infinite volume-limit, one obtains the IR exponent ac predicted by analytic 
studies [|,|,|]. 

From the above results, it seems necessary to extend present numerical simulations to very 
large lattice volumes in order to verify if the agreement obtained between lattice data and analytic 
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results for the two-dimensional case [13] applies to the 3d and 4d cases as well. Here we present 
extensive simulations in three and in four dimensions, for the SU (2) case, using huge lattices. In 
particular, we considered lattice sides = 140, 200, 240 and 320 in 3d at j8 = 3.0 and = 48, 56, 
64, 80, 96 and 128 in 4d at j8 = 2.2. In 3d, the number of configurations was about 630, 525, 350 
and 45, respectively for the four lattice sizes, both for the gluon and for the ghost propagators. In 
the 4d case we have considered for the gluon propagator 168 configurations for V = 128^ and about 
250 configurations for the other lattice sizes. For the ghost propagator we have 2 1 configurations 
for the largest volume and about 100 in the other cases. For the inversion of the Faddeev-Popov 
matrix we used the so-called point-source method [^, ^]. Note that the lattice volumes 320^ at 
j8 = 3.0 and 128"^ at jS = 2.2 correspond, respectively, ioV ^ (85fm)3 and V w {llimf. (See 



Refs. [14, 28] for details about how the physical lattice spacing a has been set in the two cases.) 
Also note that all our runs are in the scaling region. 

These simulations have been done in the IBM supercomputer at USP. This machine has 112 
blades with 2 dual-core PowerPC 970 2.5GHz CPU's, a Myrinet network and about 4.5 Tiiops 
peak-performance (occupying position number 363 in the TOP 50 list of November 2006). 

2. Results 

Considering the new data produced in the 3d case together with old data from Refs. [ |T4| , p^ ], 
we tried an extrapolation to infinite volume for the gluon propagator D{0) as a function of the 
inverse lattice side 1/L. The data for the propagator have been renormalized following Ref. []T3|]. 
Results are shown in Fig. |l|. It is clear from the plot that there is no sign of a propagator going 
to zero as L goes to infinity. Moreover, the data show a behavior of the type D(0) ~ 1/L and an 
infinite-volume extrapolation given by D{0) w 2 GeV^^. This implies c?/) = 1. 

In the 4d case (see Fig. ^), even considering very large lattice volumes and relatively large 
statistics, one cannot see a clear sign of a gluon propagator D{p^) decreasing at small momenta. 



Similar results have been presented at this conference by other groups [ |25| , |29| ]. Clearly, also in 
this case the data suggest ajj = I. On the other hand, violation of reflection positivity is confirmed 
and the gluon propagator, considered as a function of the spatial separation s, becomes negative at 



s « Ifm, in agreement with Ref. [3C]. 

We have also tried to estimate the IR exponent ac for the ghost propagator (in the 3d and 4d 
cases) using the Ansatz G{p) = cj /j^^^+'^g) ^nd considering, for each lattice volume, either the two 
smallest nonzero momenta or the third and fourth smallest nonzero momenta. Results are reported 
in Tables |l|. As one can see, this IR exponent seems to go to zero as the infinite-volume limit 



is approached, in agreement with []16||. One should however notice that, for p 500 MeV, the 



exponent qq is about 0.3, also in agreement with Ref. []16[]. 
3. Conclusions 

The above results leave us with several open questions. From the lattice point of view, one 
should of course investigate if Gribov-copy effects and/or finite-volume effects could explain our 
results. Let us recall that an improved gauge-fixing method [^] seems capable of reducing finite- 
volume effects for the gluon propagator by enlarging the set of allowed gauge transformations. At 
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Figure 1: Renormalized gluon propagator at zero momentum /3a^D(0) (in GeV ^) as a function of the 
inverse lattice side 1/L (in GeV) and extrapolation to infinite volume. The fit is given by b + c/U with 
e = 1.04(5) and = 2.05(5) GeV-^. 





Figure 2: Unrenormalized gluon propagator a^D{p^) (in GeV ^) as a function of the momentum p/a (in 
GeV) for lattice volumes V = 80"^ (left) and V = 128"* (right) at j3 = 2.2. 
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140^ 0.073(4) 

200^ 0.051(3) 

240^ 0.003(3) 

320^ -0.021(9) 



A^ 



Ug 



48^ 0.093(7) 

56^ 0.063(6) 

64^ 0.049(9) 

80^ 0.052(5) 

112^ 0.038(6) 

128^ 0.016(5) 



ClG 



140^ 0.13(2) 

200^ 0.06(2) 

240^ 0.10(2) 

320^ 0.01(5) 



Ug 



0.19(4) 
0.18(4) 



56^ 

64"^ 0.17(4) 

80^ 0.15(2) 

112^ 0.10(7) 

128^ 0.06(3) 



Table 1: Table for the ghost propagator IR exponent ac, in the 3d and 4d cases, obtained using either the 
two smallest nonzero momenta (left) or the third and fourth smallest nonzero momenta (right). 



the same time, one needs to reconcile the above results with the non-renormalizability of the ghost- 
gluon vertex [32] and with the suppression of D{p^) when considering simulations in the strong- 
coupling regime [p3|], in the interpolating gauge (or A -gauge) |Q and in Coulomb gauge [^]. 
From the analytic point of view, it may seem necessary to reconsider partially the conventional 
confinement scenarios [|l|, ^, ^, ^ ^ ^, ^ |8[ ^ 10, |n|, 12] discussed in the Introduction. One 
should of course recall that there are different solutions of Dyson-Schwinger equations for gluons 
and ghosts in Landau gauge [36, |^]. In particular, the results obtained in Ref. [36] support a 
finite non-zero gluon propagator and an essentially tree-level ghost propagator at small momenta. 
Similar results are obtained in Ref. [38]. Phenomenological tests [ ^ ] also seem to favor D{0) > 0. 

We believe that a clarification of the present status of the Kugo-Ojima/Gribov-Zwanziger sce- 
nario will probably require new ideas and new methods, both for analytic and numerical studies, 
and that a key point will be a better understanding of the gauge interpolating between the Landau 
and the Coulomb gauge [pO|]. 
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